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ABSTRACT 

The  recent  development  of  channelled  spectropolarimetry  presents  opportunities  for  spectropolarimetric  mea¬ 
surements  of  dynamic  phenomena  in  a  very  compact  instrument.  We  present  measurements  of  stress-induced 
birefringence  in  a  simple  plastic  by  both  a  reference  rotating-compensator  spectropolarimeter  and  a  channelled 
spectropolarimeter.  The  agreement  between  the  two  instruments  shows  the  promise  of  the  channelled  tech¬ 
nique  and  provides  a  proof-of-principle  that  the  method  can  be  used  for  a  very  simple  conversion  of  imaging 
spectrometers  into  imaging  spectropolarimeters. 
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1.  INTRODUCTION 

The  need  to  measure  retardance  accurately  is  common  in  optical  systems,  and  the  many  different  techniques 
employed  each  offer  their  own  unique  strengths  and  weaknesses,  with  the  requirements  of  a  given  application 
informing  the  suitability  of  the  polarimeter  design.  For  reliability  and  accuracy,  a  common  choice  is  to  use  a 
“time-sequential”  instrument  with  rotating  polarization  elements.1  For  high-speed  polarimetry,  electro-optic  or 
photoelastic  modulators  are  often  used  in  place  of  a  rotating  element.2  For  instantaneous  acquisition  of  polar¬ 
ization  information,  division-of-aperture3  or  division-of-amplitude4  designs  are  often  employed.  The  technique 
of  channelled  spectropolarimetry,6  while  not  widely  known,  offers  advantages  of  its  own  and  is  well  suited  to 
measurements  of  the  wavelength-dependence  of  birefringence.  Photoelasticity  applications  can  benefit  substan¬ 
tially  from  obtaining  spectral  information,6  and  we  present  a  method  which  uses  the  derived  birefringence  to 
avoid  some  of  the  ambiguities  of  photoelastic  fringe  analysis. 

2.  CHANNELLED  SPECTROPOLARIMETRY 

The  hardware  configuration  of  a  channelled  spectropolarimeter  (CHSP),  illustrated  below  in  Fig/^cbmprises 
two  thick  retarders,  Ri  and  R2,  oriented  respectively  at  0°  and  45°  to  horizontal,  followed  by  an  analyser  A  and 
the  spectrometer.  The  electric  field  of  light  passing  through  the  system  can  be  most  conveniently  represented  in 
terms  of  the  Stokes  vector  representation 
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where  each  of  the  three  polarization  components  (sj,  s2,  S3)  are  intensities  given  by  differences  in  orthogonal 
elements  of  the  beam.  Using  Mueller  calculus,  we  can  then  calculate  the  spectral  polarization  state  of  the  beam 
by  cascading  Mueller  matrices  for  each  optical  element.  Using  to  represent  a  linear  analyzer  oriented  at 

Further  author  information:  (Send  correspondence  to  E.L.D.)  E-mail:  eustace@u.arizona.edu 


20060824226 


1 


Distribution  Statement  A  -  Approved  for 
public  release,  distribution  is  unlimited. 


t 


Polarizer 

Tungsten-  |  QWP 
Halogen 
Lamp 


V 


R,  R2  A 


L 

Ocean  Optics 
USB2000 

Spectrometer 

PSG 


Sample  L 


CHSP 


Figure  1.  A  diagram  of  the  experimental  layout,  using  a  broad  spectrum  polarization  state  generator  (PSG),  the  sample, 
and  a  channelled  spectropolarimeter  (CHSP). 


angle  <j>  to  the  horizontal,  and  R(<f>,  S)  for  a  linear  retarder  oriented  at  angle  <f>  with  retardance  6,  then  the  Mueller- 
calculus  of  the  system  will  take  an  input  polarization  state  sin  and  produce  an  output  (before  the  spectrometer) 
of 


Bout  =  A{ 0)  •  fl(f ,  S2)  ■  R( 0,  Ji)  •  Bin  =  1/2 


/ sq  -(-  Si  cos (S2)  +  32  sin(5i)sin(<52)  -  s3  cos (5X)  sin(62)'' 
so  +  si  cos(<52)  +  32  sin(5i)  sin(62)  —  S3  cos(5i)  sin(<52) 

0 

V  0  / 


(2) 


The  detector  itself  is  directly  sensitive  only  to  the  total  intensity-component  (s0)  of  the  Stokes  vector,  and  thus 
acts  as  a  device  with  a  Mueller  matrix  of  (1  0  0  0).  Multiplying  the  above  output  Stokes  vector  with  the  detector 
sensitivity  matrix  gives  an  output  of  the  form 


2 P(cr)  =  s0  +  Si  coh(<52)  +  s2  sin^^sin^)  -  s3  cos(^i)  sin(<52)  . 


(3) 


From  this  result,  we  can  see  that  the  four  Stokes  components  are  each  encoded  as  sinusoidal  modulations  in  the 
intensity  spectrum  reaching  the  detector;  and,  since  the  retardances  61  and  S2  are  both  wavenumber-dependent, 
the  encoded  power  spectrum  itself  is  wavenumber-dependent  (see  Figs  2&3  for  examples  of  the  unencoded  and 
encoded  spectra).  Rewriting  the  above  equation  in  terms  of  complex  exponentials,  we  can  gather  terms  with  like 
carrier  terms,  giving  the  result: 

P{<t)  =  5*0+  \si{eiSl  +e~iSi) 

+|[(s2  -  is3)e«s*+Sl'>  +  (s2  +  is3)e-i^+Sl'>  +  (-s2  -  ts3)c<(is_il)  +  (-*2  +  is3)e-'<^-^)]  ,  (4) 

recalling  that  51,53,53,6!,  and  S2  are  all  functions  of  wavenumber  cr.  The  seven  terms  in  the  equation  above 
represent  the  seven  “channels”  into  which  the  polarization  information  is  encoded.  If  we  now  take  the  autocor¬ 
relation  function  (Inverse  Fourier  Transform)  of  this  power  spectrum,  we  can  then  see  the  resulting  7-channel 
distribution,  shown  in  Figs  4&5. 

Before  discussing  how  to  reconstruct  the  Stokes  spectra  from  the  autocorrelation  function,  it  is  first  useful 
to  have  a  concrete  idea  of  the  meaning  of  the  function’s  abscissa.  If  we  first  consider  the  case  of  input  light 
polarized  along  the  horizontal  axis,  then  we  can  note  that  this  light  will  be  unaffected  (other  than  by  an  absolute 
phase)  by  J?i,  whose  fast  axis  is  aligned  to  the  horizontal.  As  this  light  passes  into  the  R2,  however,  it  will  be 
separated  equally  into  two  orthogonal  elements,  one  aligned  along  R3' s  fast  axis  and  the  other  along  its  slow 
axis,  such  that  the  two  components  of  the  originally  horizontally-polarized  beam  will  experience  a  relative  path 
difference  equal  to  J?2’s  retardance.  As  this  light  passes  through  the  analyzer  it  will  experience  an  attenuation 
dependent  upon  the  polarization  angle,  and  since  the  retardance  is  heavily  wavelength-dependent  in  a  thick 
retarder,  this  polarization  angle  will  vary  rapidly  with  wavenumber,  creating  a  sinusoidally-varying  attenuation 
of  the  transmitted  beam  with  respect  to  wavenumber.  Thus,  the  abscissa  of  the  autocorrelation  function  readily 
appears  as  the  optical  path  difference  experienced  by  the  two  orthogonal  elements  comprising  each  individual 
Stokes  component  (i.e.Io  and  /90  in  the  case  of  si). 
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Figure  2.  The  measured  channelled  spectrum  for  an  un¬ 
polarized  input. 


Figure  4.  The  magnitude  of  the  channelled  spectrum 
autocorrelation  function  for  a  complex  input  polarization 
state. 


Figure  3.  The  measured  channelled  spectrum  for  a  com¬ 
plex  input  polarization  state. 
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Figure  5.  Magnification  of  the  central  region  of  Fig.  5, 
with  Fourier  windows  superimposed  (and  rescaled  for  vis¬ 
ibility). 


The  prototype  instrument  for  these  measurements  uses  a  pair  of  retarders  whose  thickness  ratio  is  3:1, 
producing  the  distribution  of  bands  in  the  autocorrelation  function  as  shown  in  Figi  5.  It  is  known  that  a 
thickness  ratio  of  1:2  actually  provides  a  more  efficient  use  of  the  spectral  bandwidth,7  lacking  the  wasted  region 
between  bands,  yet  the  configuration  in  use  here  is  still  functional.  Additionally,  the  spectrometer  chosen  for  the 
instrument  provides  a  wavelength  resolution  which  is  considerably  finer  than  the  modulations  introduced  into 
the  encoded  spectrum,  resulting  in  some  wasted  bandwidth  at  the  outer  portions  of  the  autocorrelation  function 
which  are  apparent  in  Fig.  4.  An  upgrade  of  the  prototype  will  involve  use  of  thicker  retarders  with  the  ideal 
thickness  ratio. 

With  each  of  the  Stokes  components  separated  in  the  OPD-domain,  we  can  then  mask  off  a  given  component 
with  a  Fourier  windowing  function  and  inverse  transform  to  derive  its  distribution  in  the  spectral  domain  (a 
“Stokes  spectrum”).  A  result  of  the  use  of  Fast  Fourier  Transform  routines  in  the  reconstruction  is  the  creation 
of  artifacts  at  both  the  low- A  and  high-A  ends  of  the  spectral  region.  By  proper  selection  of  the  Fourier  windowing 
function,  these  artifacts  can  be  restricted  to  the  edge  region  such  that  truncation  will  result  in  a  minimal  loss  of 
data. 
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3.  ROTATING-COMPENSATOR  SPECTROPOLARIMETRY 

As  a  reference  against  which  to  compare  the  prototype  CHSP  instrument,  we  constructed  a  standard  polarimeter 
design  which  we  refer  to  as  a  rotating-compensator  spectropolarimeter  (RCSP),  consisting  of  a  rotating  quarter- 
wave  plate  (oriented  at  various  angles  <j>),  followed  by  an  analyzer  and  spectrometer.  Mueller  calculus  for  this 
kind  of  system  gives  the  form  of  light  reaching  the  detector  as 


Bout  ~  A(0)  •  R(tf>,  5)  -  Sin 


=  1/2 


^  Sq  +  si  [cos2  (2<j>)  +  sin2(2</>)  cos  5]  +  s2[sin(2<£)  cos(2</)(l  —  cos  5)]  —  s3[sin(2</>)  sin 
s0  +  si[cos2(2<^)  +  sin2  (2<f>)  cos  5]  +  S2[sin(2ip)  cos(2^)(l  —cos  5)]  —  s3  [sin(2<£)  sin  5] 

0 

V  0  / 


(5) 


The  detector  output  will  thus  take  the  form 

2 P(u)  =  s0  +  si  [cos2(2 <j>)  +  sin2(2i j>)  cos  5]  +  s3{(l  -  cos  $)sin(2<j>)  cos(2<£)]  —  «3 [sin(2</>)  sin  <5]  .  (6) 

The  retarder  used,  while  nominally  a  quarter- wave  achromat,  has  significant  dispersion  of  retardation,  and  so 
<5  is  a  function  of  wavenumber.  The  orientation  of  the  fast  axis  <f>  also  has  a  slight  but  significant  wavelength 
dependence,  so  that  both  of  the  two  retarder  parameters  had  to  be  carefully  calibrated  to  achieve  accurate 
results.  Rewriting  the  functional  form  of  P  in  terms  of  the  orientation  angles  used  and  the  intensities  measured 
gives  P(<j>)  —  W(s,  <f)  ■  8in,  where  W  is  the  “system  matrix”  connecting  measured  intensities  with  the  input 
Stokes  components.  Written  out  in  full,  this  equation  has  the  form 


(W{« o,<t>o)  W(sU(M  WM o)  W(s3,4> o)\ 

P{4>  i) 

W(.s0,<6i)  Wiauk)  W(S2,<h) 

P(h) 

= 

W(s0,<h)  W(su<h)  W(s2,<h)  WM 2) 

P(4*) 

W{s0,<j>3)  W{auti)  W(s2,<h)  W{a3,<h) 

\P(MJ 

\WMt)  W(su<i>4)  W(s2,<t>i) 

Si 

«2 

\  3/ 

(7) 


In  this  form  it  is  easy  to  see  that  the  system  matrix  can  be  pseudo-inverted  upon  insertion  of  the  known 
orientation  angles  and  retardance.8  Obtaining  the  pseudoinverse  W+,  we  can  multiply  W+  into  the  above 
matrix  equation  giving  sin  —  W+(s,<j>)  ■  P(4>),  our  desired  result.  Performing  this  analysis  at  each  wavelength 
allows  the  reconstruction  of  the  full  Stokes  spectra.  Although  the  RCSP  analysis  can  be  performed  using 
Fourier  techniques,  we  find  the  pseudoinverse  method  more  straightforward  when  dealing  with  the  wavelength- 
dependence  of  5  and  <j>  in  the  waveplate. 

Having  our  prototype  and  reference  instruments,  we  can  use  the  polarization  state  generator  to  send  in  known 
polarization  states,  and  compare  the  fidelity  of  SOP  reconstruction.  Fig.  6  shows  two  such  comparisons,  from 
which  several  features  are  noticeable:  (1)  the  edge  artifacts  in  the  CHSP  reconstruction  can  be  entirely  eliminated 
by  truncating  the  spectral  region  of  interest  from  400-700 nm  to  415-685  nm.  (2)  While  the  reconstructions  of 
90°  LP  show  excellent  agreement,  there  is  a  substantial  discrepancy  in  that  for  RCP,  which  we  suspect  to  be  due 
to  digitization  effects  and  aliasing  resulting  from  the  use  of  only  a  minor  portion  of  the  total  available  spectral 
bandwidth. 


4.  MEASURING  RETARDANCE 

To  evaluate  the  performance  of  each  polarimeter,  we  can  now  generate  a  spectrum  with  a  known  SOP,  pass 
it  through  a  sample,  and  measure  the  output  SOP.  To  perform  a  linear  retardance  measurement  for  the  two 
polarimeter  instruments  we  present  here,  it  is  ideal  to  use  a  broad-spectrum  source  which  outputs  circularly- 
polarized  collimated  light.  Into  this  beam  we  place  the  birefringent  material  at  normal  incidence.  The  spec¬ 
tropolarimeter  is  then  placed  to  measure  the  output  SOP  across  the  entire  wavelength  region  of  interest  (in  this 
case  415-685  nm). 

In  the  Poincar4  Sphere  representation  of  the  Stokes  vector,  the  action  of  a  linear  retarder  on  the  input  Stokes 
vector  is  to  rotate  it  by  an  angle  6  (the  magnitude  of  the  retardation  in  radians)  where  the  axis  of  rotation  is 
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Figure  6.  RCSP  and  CHSP  reconstructions  for  (a)  horizontally-polarized  input,  and  (b)  RCP  input  (si  =  solid,  e2  = 
dashed,  s3  =  dotted).  The  edge  artifacts  appear  only  in  the  CHSP  data. 

twice  the  orientation  angle  (j>  of  the  waveplate  fast  axis  with  respect  to  horizontal.  Since  the  retardation  of  a 
birefringent  material  is  in  general  wavelength-dependent,  the  result  will  be  a  “Stokes  snake”  strung  across  the 
surface  of  the  Poincard  Sphere.  (For  an  illustration  of  this,  see  Fig.  7.)  So  far,  our  analysis  does  not  yet  include 
the  effect  of  multiple-reflections  inside  the  sample,  which  can  cause  a  modulation  in  the  retardation  curve  with 
respect  to  wavelength.9 


Figure  7.  Representation  on  the  Poincare  Sphere,  showing  the  locus  of  points  on  the  great  circle  generated  by  a  retarder 
oriented  at  angle  <j>  on  an  RCP  input. 

The  projection  of  the  Stokes  snake  onto  the  Poincar4  Equator  (the  «i-s2  plane)  produces  a  linear  distribution 
of  Stokes  vectors  at  right  angles  to  the  azimuthal  (a  =  2$  axis  of  the  sample  (shown  figuratively  in  Fig.  7; 
actual  measurements  are  shown  in  Fig.  8).  A  linear  fit  to  this  distribution  allows  us  to  establish  the  azimuthal 
angle  a  with  respect  to  the  s2  axis.  Half  this  angle  is  then  the  orientation  angle  <j>  of  the  fast  axis  of  the  retarder. 

As  the  Stokes  snake  wraps  around  the  surface  of  the  Sphere,  it  traces  out  a  circle.  If  the  input  polarization 
vector  is  circularly  polarized  (i.e.  at  one  of  the  poles  of  the  Sphere),  then  the  trace  is  a  great  circle  of  the  Sphere, 
which  maximizes  the  information  content.  Knowing  the  azimuth  angle  from  the  step  above,  we  can  rotate  the 
Poincard  Sphere  about  the  s3-axis  such  that  the  new  coordinate  system  (s^,  s'2,  s3)  has  its  s'raxis  coincident  with 
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Figure  8.  The  “Stokes  snake”  projection  onto  the 
Poincarfi  Equator.  The  degree  of  deviation  of  the  data 
shown  from  linear  is  more  likely  a  property  of  the  complex 
photoelastic  response  of  the  sample  than  a  reconstruction 
artifact,  as  it  shows  up  clearly  in  both  instruments. 


Figure  9.  The  circle  on  the  Poincare  Sphere  traced  out  by 
the  input  RCP  Stokes  spectrum  passing  through  a  linear 
retarder.  s'2  is  the  s2  coordinate  axis  rotated  through  the 
azimuth  angle  a. 


the  azimuthal  angle  a.  In  this  new  coordinate  system,  we  can  then  plot  the  data  points  on  the  (sj,  S3)  plane, 
producing  an  arc  which  may  be  part  of  a  circle,  or  perhaps  a  circle  which  overlaps  itself  (see  Fig.  9).  We  can 
use  the  relation  tan(2S)  —  s'2/.%  to  determine  the  angle  of  each  given  vector  and,  knowing  the  input  SOP,  we 
have  the  relative  retardation  5,  folded  by  modulo  2n.  In  order  to  determine  the  absolute  retardation,  we  need  to 
supply  a  missing  constant  n  number  of  rotations, 

The  relation  of  retardance  to  wavenumber  follows  the  relation  6(0)  —  2nd  Ana,  where  d  is  the  physical 
thickness  of  the  retarder  and  An  is  its  birefringence.  Here  we  have  a  linear  equation  in  sigma  implying  that 
5(a)  =  0  when  a  —  0,  which  provides  the  missing  constant  needed.  This  only  holds  true  when  An  is  a  constant, 
but,  for  materials  which  are  transparent  in  the  visible  region,  the  assumption  of  a  constant  is  not  a  bad  one  and 
the  resulting  model  can  produce  good  results.  Having  our  retardance  values,  we  can  then  measure  the  physical 
thickness  of  our  sample  to  derive  the  birefringence,  An  =  X6/2nd 


Figure  10.  The  reconstructed  absolute  retardance  values  for  the  sample. 

For  monochromatic  polarimeters,  multiple  measurements  are  normally  necessary  to  unambiguously  determine 
the  azimuthal  axis  of  the  sample  under  scrutiny,  but  gathering  the  polarization  data  over  a  wide  spectrum 
eliminates  the  ambiguity.  For  the  RCSP  design,  it  is  of  cour&e"St$U  necessary  to  acquire  multiple  measurements 
to  obtain  the  Stokes  snake,  but  the  CHSP  is  able  to  acquiWthisyinformation  in  a  single  snapshot.  With  the 
CHSP’s  acquisition  time  limited  by  the  speed  of  the  spectrofrtefer,  this  technique  is  not  necessarily  fast,  but 
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it  does  allow  the  instrument  to  capture  dynamic  phenomena  on  a  substantially  shorter  timescale.  Thus,  while 
slower  than  many  monochromatic  techniques,  the  channelled  spectropolarimeter’s  measurement’s  access  to  the 
orientation  angle  can  provide  a  distinct  advantage  for  some  applications. 

5.  PHOTOELASTIC  MODEL 

For  many  isotropic  materials,  stresses  applied  to  the  material  along  a  single  axis  alter  it  into  a  uniaxial  one,  where 
the  induced  birefringence  axis  coincides  with  that  of  the  applied  stress.  Thus,  photoelastic  experiments  are  an 
instance  of  retardation  measurement.  In  the  real  world,  no  material  is  completely  isotropic,  and  its  response  to 
applied  stress  can  be  much  more  complex  than  assumed  in  this  model;  whereas  some  materials  such  as  epoxy 
resin  can  use  the  simplified  model  to  a  high  level  of  approximation,  in  the  experiment  we  present  here,  the  sample 
material  is  acrylic,  whose  photoelastic  response  fits  the  simplified  model  only  with  difficulty. 

Acrylic  in  general  contains  large  residual  stresses  frozen  into  the  material  during  the  manufacturing  process. 
Thus,  its  photoelastic  response  to  an  applied  stress  will  be  a  linear  superposition  of  both  the  frozen-stress  and 
applied-stress  terms.10  In  addition,  its  plasticity  gives  the  material  a  “memory”,  such  that  after  applying  stress 
it  may  require  a  substantial  relaxation  time  to  return  to  its  original  state.  Even  after  relaxation  the  material 
might  not  return  to  its  original  state,  thus  every  sequence  of  measurements  must  be  preceded  by  a  measurement 
of  the  material  without  applied  stress  in  order  to  obtain  a  baseline.  Additionally,  the  plasticity  of  the  sample 
induces  a  time-dependent  photoelastic  response:  a  continuously  applied  constant  stress  will  induce  a  steadily 
changing  birefringence  in  the  sample.  Thus,  the  measurements  carried  out  here  we  performed  as  quickly  as 
possible  to  minimize  the  time-dependence  of  the  material. 

During  the  experiment,  we  applied  steadily  increasing  linear  stresses  to  the  sample,  taking  data  at  each  step. 
Each  Stokes  spectrum  obtained  thus  represents  the  polarization  state  of  the  beam  exiting  from  the  sample, 
averaged  across  the  beam  diameter.  The  current  experiment  therefore  assumes  that  the  material  properties 
are  uniform  across  the  beam  aperture.  Future  research  will  be  directed  at  implementing  the  channelled  spec- 
tropolarimetry  technique  to  an  imaging  system  which  will  obviate  the  need  for  this  assumption.  If  we  treat  the 
sample  as  a  linear  retarder  of  unknown  retardation  amplitude  and  unknown  orientation,  we  can  then  use  the 
“Photoelastic  Law”  to  derive  the  differential  stress  experienced  by  the  sample11: 

As  =  An/C  ,  (8) 

where  As  is  the  differential  stress  (expressed  in  units  of  force  per  unit  area)  and  C  is  the  stress-optic  coefficient 
for  the  sample  material.  Thus,  from  the  birefringence  of  the  sample  derived  in  the  previous  section,  we  can 
either  use  the  known  applied  stresses  to  determine  the  constant  C  of  the  sample  material,  or,  if  the  constant  C 
is  already  known,  use  the  result  as  a  direct  measure  of  the  accuracy  of  the  procedure.  In  the  case  of  plastics,  the 
optical  constant  varies  greatly  from  sample  to  sample  (even  using  the  same  nominal  chemical  composition)  even 
if  its  time  dependency  is  ignored.  The  results  of  our  experiment,  for  both  the  RCSP  and  CHSP  results  using  the 
same  sample,  are  shown  in  Fig.ll. 

Because  of  the  large  residual  stresses  frozen  into  our  sample,  it  is  natural  to  expect  that  the  measurement 
of  stress  orientation  angle  <j>  will  not  be  a  constant  with  applied  stress.  Only  for  large  values  of  applied  stress 
does  the  residual  stress  become  ignorable,  but  large  stresses  are  problematic  for  the  simple  model  applied  here. 
Nevertheless,  for  a  stress  applied  at  roughly  horizontal  to  the  optical  system  (0°),  our  fitted  results  showed  a 
stress  orientation  angle  of  3.6°for  the  RCSP  and  8.7°for  the  CHSP  reconstructions.  The  discrepancy  between 
the  two  measurements  of  <f>  is  directly  linked  to  the  artifacts  (indicated  in  Fig.6)  in  the  CHSP  reconstruction  for 
circularly  polarized  light. 

The  data  presented  here  is  only  for  very  low  stresses,  in  which  the  photoelastic  response  of  the  material  can 
be  approximated  with  a  constant  stress-optic  coefficient  C;  with  larger  stress,  the  coefficient  begins  to  decrease, 
which  in  the  birefringence  data  is  seen  as  a  levelling-off  towards  a  constant  value.  Given  the  above  measurements 
of  photoelastic  response,  we  calculate  a  stress-optic  coefficient  for  the  sample  of  1.13  x  10~12cm2/dyne. 
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Figure  11.  The  sample  birefringence  at  X  =  550  nm  as  a  function  of  applied  stress. 

6.  CONCLUSION 

The  advantages  of  the  channelled  spectropolarimetry  technique  include  its  simplicity  of  construction  and  its 
minimal  use  of  moving  components,  which  can  make  it  very  robust  and  improve  its  acquisition  speed.  The 
degree  of  agreement  between  the  RCSP  and  CHSP  instruments  shown  here  indicates  the  promise  of  the  channelled 
spectropolarimetry  technique,  despite  the  considerable  room  for  improvement  in  our  current  instrumentation  and 
reconstruction  methodology. 

Our  lab  has  for  some  time  now  been  working  on  developing  novel  imaging  spectrometer  designs,  especially 
for  applications  which  require  the  measurement  of  dynamic  phenomena.  Channelled  spectropolarimetry  presents 
obvious  advantages  for  the  extension  of  such  a  system  to  polarimetry  without  some  of  the  difficult  issues  expe¬ 
rienced  by  other  techniques  (one  example  of  which  would  be  the  co-registration  of  images  in  a  multiple-detector 
system).  By  extension  of  the  photoelasticity  experiments  shown  here  to  an  imaging  system,  we  hope  to  mea¬ 
sure  dynamic  phenomena  such  as  the  propagation  of  cracks  or  response  to  impacts  —  measurements  which  can 
profit  from  the  instrument’s  ability  to  simultaneously  map  the  spatial  distribution  of  both  stress  magnitude  and 
orientation. 
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